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Abstract 



> 

| We calculate B^ and to -^(1430) transition form factors, and study 

semileptonic decays of Br s \ and Dr s \ — > Kq(1430)£u based on QCD sum rule. 
Measuring these semileptonic decays with high statistics will give valuable 
, information on the nature of light scalar mesons. 

o 

^ . PACS numbers: 12.38.Lg, 13.20.Fc,13.20.He 

Semileptonic decays of B and D mesons are important for studying quark flavor mixing 
and extracting the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, because strong 
interactions involved in these processes are simpler than that of hadronic decays. The 
^ ■ strong binding effects can be parameterized as transition form factors, which can be cal- 

culated by Lattice QCD P, QCD sum rule [2J [3 |1J |S] , light-cone sum rule [H], or by quark 
model [Zj , light-front approach [B] , and recently by large-energy and heavy-quark-effective 
theory 0. Experimental data on semileptonic decay of B and D mesons can be used 
to test these theoretical method of treating nonperturbative dynamics. Transition form 
factors inducing semileptonic decays not only depend on dynamics of strong interactions 
between quarks in the initial and final hadrons, but also on the structure of the hadrons 
involved in the semileptonic decays. 

The structures of scalar mesons are long-standing problems in particle physics. A large 
amount of scalar mesons have been found in experiment [TJ3- They include a [or /o(600)], 
/o(980), / (1370), / (1500), / (1710), a (980), a (1450), k, #£(1430), etc. [TT]. At least 
they can be divided into two flavor nonets, one below or near lGeV, the other above 
lGeV. The structure of scalar mesons is still not well established theoretically ^21- bi 
the literature, many suggestions are discussed such as qq, qqqq and meson-meson bound 
states. Among the controversy, almost every model on scalar states agrees that #^(1430) 
is dominated by su or sd state. Recent analysis from QCD sum rule agrees with this 
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picture ^3]. The result of QCD sum rule favors that Kq(1430) is mainly sq bound state. 
The predicted mass of sq + scalar bound state is consistent with the mass of 1^(1430). 

To investigate the structure of scalar meson, a large amount of experimental data 
and theoretical studies are necessary, scalar mesons are produced in ttN scattering, pp 
annihilation, decays of heavy flavor mesons, etc.. Recently, hadronic D decays involving 
scalar mesons are studied in generalized factorization approach ^1] . Hadronic decay of D 
meson is complicated, it suffers from final-state interactions which are difficult to control 
theoretically. Compared with hadronic decays, semleptonic decays are more simpler. They 
are free from final-state interactions. All nonperturbative binding effects are parameterized 
in the transition form factors. 

In this work, we shall study semileptonic decays of B( s -> and Dr s -> meson involving 
1£q(1430). In Ref.|15j. semileptonic decays of D meson into a, k are studied in QCD sum 
rule, where the scalars a and k are treated as qq bound states, which is different from 
the picture that scalars below 1 GeV are predominantly multiquark states. According 
to our previous study in QCD sum rule, 1^(1430) is dominantly the ground state of 
sq scalar channel. We treat Kq(1430) as sq bound state in this work. The transition 
form factors of Br s -\ and Dr s \ — > .^(1430) are calculated in QCD sum rule. Then they 
are used to study semileptonic B and D decays. The semileptonic decay modes of B^ 
and D {s} mesons involving i^(1430) include 5° -> K^(U30)~£u, D° K$(U30)-£v, 
D + -» K^(U30)°£u, and D+ -> i^ *(1430) ^, where i is the lepton e or \i. Compared to 
the large mass of B and D meson, the mass of the lepton e or [i is dropped. 

I I 
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Figure 1: Semileptonic decays of and D( s ) involving Kg (1430). (a) for B® — > 
K*(U30)-£u, (b)_for D° -> K*(U30)~Iu, (c) for D + -> K*(U30)°£v, and (d) for 
D+ -> ^(1430)°^ 

The Feynman diagrams for the semileptonic B/ s ) and Dr s ^ decays into -fCg(1430)£^ 
are shown in Fig-ljljl. In the limit of flavor SU(3) symmetry, the amplitudes of the three 
decay modes D° -» K$(U30)-£v, D+ -> K^(U30)°£iy and Df -» K$(U30)°£v should be 
the same. However, flavor SU(3) symmetry is not an exact symmetry, because the mass 
of s quark is much larger than that of u and d quarks. To investigate SU(3) symmetry 
breaking effect, the mass of s quark m s is kept in our calculation. 
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The amplitude of B and D decaying into if q (1430)^ is 

A = ^V^uj^l - 15 )£(K* (1430)1^7^(1 - l5 )Q\H) (1) 

where Gf is Fermi constant, H is the heavy flavor meson B or D, Q is b or c quark, q' can 
be light quark u, d, s, and Vq ? ' is the relevant CKM matrix element. Strong interactions 
are contained in the hadronic matrix element (Kfi (1430) \q' 7^(1 — 75)<3|il). By analyzing 
the parity property of the hadrons and the current, we can know that the vector current 
does not contribute to the pseudoscalar-scalar hadronic matrix element 

{K*(U30)\q f ^Q\H)=0, (2) 

therefore, only the axial current contributes. From the Lorentz invariance, the hadronic 
matrix element of axial current can be decomposed into 

<*: o *(1430)|gV75Q|tf> = -i[(pi+ P2 rF + (q 2 ) + ( Pl - p 2 )"F_( g 2 )], (3) 

where the parameters F + (q 2 ) and F^(q 2 ) are the transition form factors, which only 
depend on the momentum transfer squared. The momenta p\ and P2 are the momenta of 
the initial and final state mesons, respectively, and q = pi — P2- 

Next we shall calculate the transition form factors in QCD sum rule method. In the 
limit of the lepton mass mi ~ 0, the form factor F_ does not contribute to the semileptonic 
decays of B and D mesons. So we shall not consider F- in this work. 

In the standard procedure of QCD sum rule for calculating transition form factors, all 
hadrons involved in the transition process are interpolated by currents with the appropri- 
ate quantum numbers relevant to the hadrons. The heavy flavor meson H (B or D) is 
interpolated by pseudoscalar current Jh( x ) = Q( x )il5Q'( x ), summation over the spinor 
and color indices being understood in this interpolating current. The final scalar meson 
Kq(1430) is interpolated by the scalar current Jr{ x ) = q{x)s(x), here q is up or down 
quark. Then the decay rate is obtained from the vacuum expectation value of time ordered 
product of the interpolating fields Jh(x), Jk(x) and the weak axial current J M = q'^^^Q- 
The vector part of the weak current is not presented here because its contribution to B or 
D — > Kq(1430) transition is zero. A three-point correlation function can be constructed 
from the vacuum expectation value of time ordered product of interpolating fields and 
axial weak current 

MpIpIi 2 ) = * 2 J d 4 xdSe l ^ x - tpi - y (0\T{J K (x)J^0)J H (y)}\0). (4) 

Phenomenologically, one can insert a full set of intermediate sates into the time ordered 
product and obtain the double dispersion relation, then the correlation function is ex- 
pressed as a sum of the contributions of the lowest lying and excited states 

u , 2 2 2, (0|^g / il^o(1430))(^(1430)|g / 1 7^75Q|g)(g|Qi75g 2 |0) 

+excited states, (5) 

where q[ and q' 2 can be u, d, s quarks, but they do not take the same flavor at the same 
time and there must be one to be s quark in q[ and q' 2 - By introducing (0\qs\KQ (1430)) = 
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fK*m K * and (H\Qij 5 q'\0) = m ™+ m ■ fn, where f K * and f H are decay constants, the 
correlation function is obtained as 

K»(pIpI<1 2 ) = 1 2 2 A K* (U30)\q' lllll5 Q\H} 

m Q + m q f 2 (mjj - Pi)(m 2 K , - pj) 

+excited states (6) 

Therefore the correlation function is related to the transition matrix element by eq.©. If 
the correlation function U^(p 2 ,p%, q 2 ) can be calculated in QCD reliably, then the tran- 
sition form factors can be extracted through eq.©. In QCD, the three-point correlation 
function can indeed be evaluated by operator-product-expansion (OPE) method in the 
deep Euclidean region 

p 2 «(m Q + m q ,) 2 , p 2 «(m s + m q ) 2 . (7) 

Although it is in the deep Euclidean region of p\ and p\ that the correlation function 
can be calculated, eq.© clearly shows that the value of n M (pf ,p%, q 2 ) can still be used 
to extract the transition form factors which is typically a Minkowskian quantity, be- 
cause the valuables p\ and p\ are clearly separated from the transition matrix element 
(Kq (1430)|^7^ 75QI-H) in the left hand of eq.®. Next we shall evaluate the correlation 
function in QCD. 

The time-ordered current operators in the three-point correlation function in eq.(j3J) 
can be expanded in terms of a series of local operators with increasing dimensions, 



i 2 J dVye v ^ x -^yT{J K {x)J^)J H (y)} 

= C ^I + C 3 ^ + C^G a a(3 G aal3 + C^a a(3 T a G aa ^ 

+ c^r^r'* + • • • , (8) 

where C^'s are Wilson coefficients, I is the unit operator, is the local Fermion field 
operator of light quarks, G^g is gluon strength tensor, V and V are the matrices appearing 
in the procedure of calculating the Wilson coefficients. Considering the vacuum expecta- 
tion value of the OPE of the interpolating field operator, we get the correlation function 
in terms of Wilson coefficients and condensates of local operators, 

nM,Pl<l 2 ) = i 2 Jd 4 xd 4 ye i ^ x -^y(0\T{J K (x)J fl (0)J H (y)}\0) 

= C 0fl + C 3 „<0|tf |0> + C 4fl (0\G a a3 G aal3 \0) + C Blt (0\*VafiT a G aa f ) *\0) 
+ C 6f ,(0\^T^T'^\0) + ■••, (9) 

The Wilson coefficients Ci's depend only on p\ and P2, according to the Lorentz structure 
of the correlation function, the result can be re-expressed by two parts in the "theoretical 
expression" 

n M Oi,p|,<7 2 ) = -if+(pi + P2) fJ , - if-(pi -P2)n, (10) 
The coefficients f± above collect all the contributions of perturbative and condensate terms 

f± = r± rt + /i 3) + /i 4) + /i 5) + /i 6) + • • • (ii) 
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where /± er * is the perturbative contribution of the unit operator, and , are 
contributions of condensates of dimension 3, 4, 5, 6, • • • operators in the OPE. The per- 
turbative contribution and gluon condensate contribution can be written in the from of 
dispersion integration, 



f± rt = I ' dsidar 



P± rt (s 1 ,s 2 ,q 2 } 



/± 4) = J dsids 2 -^ 



(si -Pi)(s 2 -pl) 

We can approximate the contribution of excited states as integrations over some thresholds 
Si and s 2 in the above equations, which is the assumption of quark-hadron duality [21 
H£j. Then equate the "phenomenological" and "theoretical" expressions of the correlation 
function in Eq. © and (|10|). we can get an equation for the form factors. But such 
equation may highly depend on the approximation for the contribution of excited states 
and the condensate of higher dimensional operators in OPE. To improve such equation, 
one can make Borel transformation over p\ and p\ in both sides, which can suppress 
the contributions of excited states and condensate of higher dimensional operators. The 
definition of Borel transformation to any function f(p 2 ) is 

(_ 2\n fin 

p 2 — > — OO 

-p 2 /n = M 2 

Some examples of Borel transformation is given in the following, 

R 1 _ 1 1 r -s/M 2 

\p 2 ,M*( s _ p2 }k (fc_i)!(M 2 ) fc 

B\ p 2 M 2(p 2 ) k = 0, for any k > 0. 

Equating the two expressions of the correlation function, subtracting the contribution 
of excited states, and performing Borel transformation in both variables p\ and p\, we 
finally obtain the sum rules for the form factor 

F + (,<) = t mQ + m f, MlMl^f^BU, (12) 
fK*rn K *m z H f H 

where Bf + denotes Borel transforming / + in both variables p\ and p 2 , Mi and M 2 are Borel 
parameters. After subtracting the contribution of the excited states, now the dispersion 
integration for perturbative and gluon condensate contribution should be performed under 
the threshold, 



2 («i -p 2 )(s2 -p 2 y 



(si -p?)(«2 -pIY 



In this work we calculate the Wilson coefficients in OPE up to the contribution of the 
operator of dimension 6. The diagrams we considered here are similar to our previous 
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work in Ref. 5; . We find again that the contributions of the diagrams for the gluon-gluon 
condensate in fig|2]exactly cancel in the sum in the case of pseudoscalar to scalar transition. 
Therefore, the gluon-gluon condensate does not contribute to B and D — ► #£(1430) 
transition form factors. 




(d) (e) (f) 



Figure 2: Diagrams for contributions of bi-gluon operator. 

The results of Borel transformed coefficient Bf + in ea. (|12j) are given in the Appendix. 

Since we do not take into account radiative corrections in QCD, we choose the param- 
eters at a fixed renormalization scale of about 1 GeV. The values of the condensates are 
taken as |21IT7|. 

(qq) = -(0.24 ± 0.01 GeV) 3 , (ss) = m 2 {qq), 
g^oTG^) = m§(M), a s {^) 2 = 6.0 x 10~ 5 GeV 6 , (13) 

ml = 0.8 ±0.2 GeV 2 . 

The quark masses are taken to be = 4.7 ± 0.1 GeV, m c = 1.3 ± 0.1 GeV, m s = 
150 ± 10 MeV and m u ~ ~ 0. The decay constants of B^, D^ s -j and #£(1430) are 
taken from the two-point QCD sum rule: f D = 200 ± 20 MeV, f B = 180 ± 30 MeV 
QS1 EEB1, and from the ratios f D Jf D = 1.19 ± 0.08, f B JfB = 1-16 ± 0.09 [Sj, we can 
obtain f Ds = 238 ± 29 MeV, f Bs = 209 ± 38 MeV. For the decay constant of #£(1430), 
we get fx* = 427 ± 85 MeV from two-point sum rule [T3] . 

The threshold parameters s\ and s° are also determined from the two-point sum rule. 
The parameter is for the threshold of B and D meson, they are taken to be s9g = 
35 ± 2 GeV 2 , s° D = 6 ± 1 GeV 2 |E|. s° is for the threshold of #£(1430), its value is 
s° = 4.4 ± 0.4 GeV 2 [E]. 

The Borel parameters Mi and M<i are not physical parameters. The physical result 
should not depend on them if the operator product expansion can be calculated up to 
infinite order. However, OPE has to be truncated to some finite orders in practice. There- 
fore, Borel parameters have to be selected in some "windows" to get the best stability of 
the physical results. We choose M\ and M2 in the region where 1) the contribution of the 
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excited states is effectively suppressed, which can ensure that the sum rule does not sensi- 
tively depend on the approximation for the excited states, and 2) the contribution of the 
condensates should not be too large, which can ensure that the contribution of the higher 
dimensional operators is small and the truncated OPE is effective. The contribution of the 
excited sates to the three-point correlation function is in the form e excited' 1,2 j where 
^excited denotes the mass of the excited state, and the series in OPE generally depend 
on Borel parameters in the denominator 1/Mi : 2, therefore, for effectively suppressing the 
contributions of the excited states and the higher dimensional operators in OPE, the Borel 
parameters should be neither too large, nor too small. We find the optimal stability 
with the requirements shown in Table ^ The regions of Borel parameters which satisfies 
the requirements of Table ^ are shown in FigEl in two-dimensional diagram of M\ and 
Mf . We find good stability of the form factors within these regions. 

Table 1: Requirements to select Borel Parameters M\ and Mf for each form factors 



Form Factors 


contribution 


excited states 


excited states 




of condensate 


of H channel 


of K* (1430) channel 


F+ sK( >(0) 


<42% 


< 10% 


<42% 


F+ K( '(0) 


< 30% 


< 17% 


< 24% 


F+ sK °(0) 


< 28% 


< 12% 


< 14% 




M| (GeV 2 ) 



(b) 



5 6 

Ml (GeV 2 ) 



Mf (GeV 2 ) 



6 

Ml (GeV 2 



Figure 3: Selected regions of Mf and Mf: (a) for F+' °; (b) for F+ ; (c) for F+' °. 

In general the higher the dimension of the operators, the smaller the relevant contribu- 
tions of the condensates. We find that, in the numerical analysis, the main contributions 
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to the form factors are from perturbative term and condensate of dimension-3 operator, 
the contributions of operator of dimension 6 are negligible. 
The final results for the form factors at q 2 = are 

F+ sK °(0) = 0.24 ±0.10, 

DK* 

F+ °(0) = 0.57 ±0.19, (14) 

F+ sK °(0) = 0.51 ±0.20. 

The error bars are estimated by the variation of Borel parameters, the variation of the 
threshold parameters Si 2 , the uncertainty of the condensate parameters, the variation 
of the quark masses and meson decay constants, and the possible a s corrections to the 
three-point function. The main contribution comes from the uncertainties of the threshold 
parameters and decay constants, which is about 10% ~ 20% of the central value, the other 
uncertainties are only a few percent. The a s corrections to two-point functions of Jb, /d 
and fx* are 10% to 30%, we estimate the a s corrections to three-point function shall be 
the same order. We take 30% as the error caused by the a s correction in this work. 

DK* D K* 

Numerically, F + and F + 3 are slightly different. The difference is completely 
caused by the SU(3) symmetry breaking effect. The numerical result shows that this 
effect is about 10%. 

There have been studies of heavy flavor meson to -fTg(1430) transition form factors in 

DK* 

the literature. In ^1], F + (0) is extracted from the data of the hadronic D decay in the 
generalized factorization model, F+ K ° (0) = 1.20 ±0.07. It is also calculated in covariant 

DK* 

light-front approach, F + (0) = 0.48 8 . Compared with these values, our present result 

DK* _ 

-f+ ° (0) = 0-57 ± 0.19 is consistent with the light-front approach. 

The physical region for q 2 in B or D — > A"q(1430)£i^ decay extends from to {tub,d — 
m K*) 2 , which is not large because of the large mass of Kg (1430). The q 2 dependence of the 
form factors in the physical range can be calculated directly from QCD sum rule, which 
has been discussed in detail in [4]. Within the physical range of q 2 in B or D to ^(1430) 
decays, there is no non-Landau-type singularity @f with the threshold parameters s\ and 
s 2 considered in this paper. Prom the sum rule for the form factors listed in the appendix, 
we see that the formulas are not singular at q 2 = 0, and q 2 appears in the numerator, 
therefore, if \q 2 \ is too large, OPE will fail. So the condition for the OPE be effective is to 
keep \q 2 \ in a moderate range. For the whole physical range of q 2 considered in this paper, 
within the selected window of the Borel parameters, the contribution of higher-dimensional 
operators is effectively small, therefore the OPE is effective. 

The form factors as functions of q 2 , normalized by the value at q 2 = GeV 2 , i 7 ±((7 2 )/i 7 ±(0) , s 
are shown in FigHJ Within the physical range, the behavior of the form factor F + (q 2 ) is 
compatible with the pole-model, 



F + (q 2 



F+(0) 



i-^/Sole' 



The result of the form factors from QCD sum rule can be fitted with the pole model. 
The fitted pole masses are, 

m Bs fo = 4.80 ±0.12 GeV, 



pole 

m^jO = 2.9 ± 0.3 GeV, (15) 

m DsK* = 196±012 GeV. 
pole 
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Figure 4: q 2 dependence of the form factors from QCD sum rule. The solid curve is for 
the central value, the dotted curve is for the errors allowed by the input parameters. 



Next we shall use the form factors calculated in QCD sum rule to study the differential 
and total decay rates of B and D — > Kq(1430)£u decays. The differential decay rate is 
calculated to be 



dq 2 



G 2 F \VQq'\ 2 

192vr 3 m^ 



2\2 



F + {q 2 Y[{m 2 H + m 2 -q 2 ) 



4m 2 H m% 



13/2 



(16) 



The values of the CKM matrix elements are taken in Wolfenstain parameterization, A = 
0.8, A = 0.22, p = 0.20, and r] = 0.34 M, which are relevant to V cs = 0.976, V cd = -0.22. 

The differential decay widths as a function of momentum transfer squared q 2 are 
shown in FigEI The decay rate of D+ -> ^(1430)°! v is one more order smaller than 
that of D° -> K *(1430)-^ and D + -> K$(UW)°iiy because of the CKM matrix element 
suppression of V C( i over V cs . After integration over q 2 in the whole physical region, we get 
the integrated decay widths 



ri.6 



T(B° S -> K£(U30)-£v) 
T(D°' + -» K$(U30)£v) 
F(D+ -» ET *(1430) ^) 

The branching ratio is defined by 

Br = T/F tota i 



1.1 ) 



x 10 



-17 



GeV, 



(2.9Zi%) x !0" GeV, 
(3.2l|°) x 10- 17 GeV. 



(17) 



(18) 



where V denotes the decay width of each decay mode, and T tota i the total decay width 
of B or D meson. The total decay width of one meson is related to its mean life time 
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dT(q 2 )/dq 2 dT(q 2 )/dq 2 
(xlO-^GeV- 1 ) (xlO^GeV- 1 ) 




0.25 0.3 

q 2 (GeV 2 ) 



Figure 5: Differential decay widths of B and D to Kq(1430)£v as function of q 2 . (a) For 
B° -» K$(U30)-£v; (b) The solid curve is for D° -» K$(U3Q)-£v and D+ -> K$(U30)°£v 
with the scale parameter A = 15, while the dashed curve is for —* Kq(1430)°£u with 
the scale parameter A = 16. 

r by Ttotai = Fi/t. The mean life times of B and D mesons are: tb s = 1-46 x 1CP 12 s, 
t d ± = 1.04 x 10~ 12 s, t d o = 0.41 x 10~ 12 s, and t Ds = 0.49 x 10~ 12 s 10 . We can obtain 
the branching ratios 

Br(B° -» Xg(1430) - A/) = (3.6±||) x 10~ 5 , 
fl r (D° ff*(1430)-/i/) = (1.8±i:g) x 10~ 4 , 

Br{D+ -> i^ *(1430) ^) = (4.6l|£) x 10" 4 , (19) 
Br(Df -» KJ(1430)°li/) = (2.4±?;§) x 10~ 5 . 

Although r(L>° -» ^(1430) - ^) = T(L> + ^(1430)°^) under the isospin symmetry, 
the branching ratios of these two decay modes are greatly different. This is due to the 
large difference of the total decay widths of D° and D mesons. 

The numerical result shows that the branching ratios of B and D to -K"q(1430)£z^ 
are rare. Measuring them in experiment needs high statistics, but these measurements 
will give valuable information on the nature of the light scalar mesons. Recently, Focus 
collaboration gives an upper limit for the ratio, r ^-t^n+v)^ ^ 0.64%, in analyzing 

the mass spectrum of D + — > K~ir + decay [19j . By considering the branching ratio 
of Br(D + — > K~ -k + n + v) = 4.00 ± 0.32% given in PDG ^U], one can get the branching 
ratio Br(D+ -> K$(U30)° fj,+v) < 2.8 x 10~ 4 . Compared with this upper limit, our result 
Br(D + -> ^(1430)°lz/) = (4.6±|;|) x 10" 4 is compatible with it. Direct experimental 
measurement of these semileptonic branching ratios and more precise theoretical prediction 
are highly desired. 

In summary, we have calculated the B and D to 1£q(1430) transition form factors 

sum rule, and studied the semileptonic decays of 
B° -► KS(1430)-tu, D° -> KS(1430)-£v, D+ -► K$(U30)°£v, and D+ -> K$(U30)°Ev. 
Measuring these decay modes with high statistics will give valuable information on the 
structure of light scalar mesons. 
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Appendix 

The analytical results of the coefficients Bf + for the form factors are given here. Two 
cases are classified, one is for B s and D s — > Xq(1430) transitions, the other for B and D 
transitions. 

(1) For B s and D s — > Xq(1430) transitions, the form factor in eq. (|12[) will be 

F + (q 2 ) H ^ = t mQ + m 2 \ MlMh m »^e m k /M hf^ K ^ (20) 

/ Kg m Kg m H s f H s 

where H s denotes B s or D s , and 

&f H.KS = £ fP+e rt + ^ + ^ + ^ (21) 

The perturbative contribution is 

.-0 „ o si _ s 2 



-(q 2 - si + s 2 )m^m s + mg[(g 2 + si - s 2 )s 2 + (q 2 - si + s 2 )m 2 ] 
+2q 2 [-s lS2 + {-q 2 + Sl + s 2 )m 2 }}, (22) 



where the lower integration limit of dsi is determined by the condition that all internal 
quarks in the perturbative diagram be on mass shell [20"] . 



r m Q 2 

*l = — 2 s 2 + m Q , 

rrin — Q 



and X = (si + s 2 - q 2 ) 2 - 4s 1 s 2 . 

The condensate contributions are 



4 



+M|m Q [-M 1 2 g 2 + (M 2 + MDm^m 2 }, (23) 



— Q 

®- = ~ 1!m| M| } ^mW^M 2 + 2M| + g 2 ) - (M 2 + Mf )m 2 Q ] 

+ (4M 2 + m^)[-M 2 g 2 + (M 2 + M 2 )m 2 Q ]m s }, (24) 



* (6) 47ra s (^ 2 — 

324.^^1 (rnfj - q 2 



2 

2\ 



+6q 2 )} + Mf[-38M? + g 2 (24M 2 2 + 5q 2 ) + 2M 2 (M 2 2 + l9q 2 )]m 2 Q 
-2[UMf - 2M 2 V + M 1 2 (12M| + 5q 2 )]m^ + (5M 2 - 4M 2 >^} 
+36M 1 6 M|m Q (3M 2 2 + g 2 - m^)m s - 18M 1 6 (4M 2 4 - M^m 2 Q )m 2 s }, (25) 



11 



(2) For B and D — > ETq(1430) transitions, the form factor in ea.(fl2"]l will be 

F + (q 2 ) H « K *o = 7 ^ MfMfe^^e^^B/f^o (26) 

where H q denotes B or D, and 

B/f 9 A '' = Bf + t + Bff Hq + Bf% q + Bffl ■ (27) 



The perturbative contribution is 



-3 ~ s '' : 



+ [(si - S2) 2 - g 2 (si + S2)]"iQ"i s + (g 2 - si + s 2 )mQm s 

+ Sl [2q 2 s 2 -{q 2 -s l +s 2 )m 2 s ]}, (28) 

The condensate contributions are 

" "2i|^^ S| <" ! «- m *>- (29) 



2\^, 1 A/f2/'A / r2 1 w2-\_2 , s#2^, _ r Ji,r2 / Ji,r2 , o)\#2 , „2\ 



+q z )m s + M${M{ + M^)m 2 Q m s + M(mQ[M£{M( + 2M 2 + <f 
-(Mi +M|)m 2 ]} (30) 



^ - - 324Mf M|K - q 2 f ^ {2M! [2M|(11M| - « V - MlV(6M| + g2) 
+M 4 (54M 2 2 + 23q 2 )]m Q m s - 2[23Mf + 7M 2 V - 2M 1 4 (3M 2 2 + g 2 ) 
+M 2 (22M 2 4 - 5Miq 2 )]m 3 Q m s - 2(M 4 + 3M\mI - 7M 2 4 )m|m s 
+(5M 2 - 4Ml)m 6 Q (rr4 - m 2 ) - Mfm 2 Q {-Ml[-38Mf + g 2 (24M 2 2 + 5q 2 ) 
+2M 2 (M 2 2 + 19g 2 )] + [58M X 4 + g 2 (27M 2 2 + 5q 2 ) + 6M 2 (2M 2 2 + 7g 2 )]m 2 } 
+m 4 3 {-2M 2 2 [13M 1 4 - 2M 2 V + M 2 (12M 2 2 + 5q 2 )} + [29M 4 - 4M 2 V 
+m 2 (27M| + 10g 2 )]m 2 } + M 1 4 {2M 2 2 [M 1 2 (72M 2 2 - 17g 2 ) - g 2 (M 2 2 + 6q 2 )} 
+ [-4m 2 (18M| - 19g 2 ) + q 2 (\2Ml + 13g 2 )]m 2 }}, (31) 
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